subroutine wmnij_t1ej_imnie(wmnij, iijka,& 
              t1, icore, fact)
!
use mod_size
use mod_ioff
use mod_orbit
use mod_iop 
implicit none
!
real*8, intent(in) :: fact
real*4, intent(in) :: wmnij(*)
real*4, intent(in) :: iijka(*)
!
real*4, intent(inout) :: t1(*), icore(*)
!
integer :: i, j, a, b, nsize, i0, i1, i2
integer, external :: idsymsoc
real*8, external :: dnrm2
i0 = 1
i1 = i0 + nvvoo*isd
i2 = i1 + nvvoo*isd
!this is to do wmnij = iijkl + P(ij)(iijka(m,n,i,e)*t1(e,j))
!===AAAA
!--iijka(MNIE)*t1(EJ)
call VpqrmVms_to_Vpqrs(iijka(ioi2(1)), t1(ifvo(1)), 1, 1, vrta, &
            popa, popa, 0, popa, popa, icore(i0), 1, fact, 0.d0, 0)
call iassymso(icore(i0), icore(i1), 1, popa, popa, popa, popa, 0)
call axpyso(nioi11, 1.d0, icore(i1), 1, wmnij(ioi1(1)), 1)
!===ABAB
!--iijka(MnIe)*t1(ej) - iijka(Mn,jE)*t1(EI)
! =iijka(MnIe)*t1(ej) + iijka(nM,jE)*t1(EI)
call VpqrmVms_to_Vpqrs(iijka(ioi2(4)), t1(ifvo(1)), 1, 1, vrta, &
            popb, popa, 1, popb, popa, icore(i1), 1, fact, 0.d0, 1)
call isymtrso(icore(i1), icore(i0), popb, popa, 1, 0, popb, popa, 1, 1, 1)
call isymtrso(icore(i0), icore(i1), popb, popa, 1, 1, popa, popb, 1, 0, 1)
!
call VpqrmVms_to_Vpqrs(iijka(ioi2(2)), t1(ifvo(2)), 1, 1, vrtb, &
            popa, popb, 1, popa, popb, icore(i1), 1, fact, 1.d0, 1)
!
call axpyso(nioi12, 1.d0, icore(i1), 1, wmnij(ioi1(2)), 1)
!
!===BBBB
!--iijka(mnie)*t1(ej)
call VpqrmVms_to_Vpqrs(iijka(ioi2(3)), t1(ifvo(2)), 1, 1, vrtb, &
            popb, popb, 0, popb, popb, icore(i0), 1, fact, 0.d0, 0)
call iassymso(icore(i0), icore(i1), 1, popb, popb, popb, popb, 0)
call axpyso(nioi13, 1.d0, icore(i1), 1, wmnij(ioi1(3)), 1)
return
end
